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Relativistic third-order MBPT is apphed to obtain energies of ions with two valence electrons 
in the no virtual-pair approximation (NVPA). A total of 302 third-order Goldstone diagrams are 
organized into 12 one-body and 23 two-body terms. Only third-order two-body terms and diagrams 
are presented here, owing to the fact that the one-body terms are identical to the previously studied 
third-order terms in monovalent ions. Dominant classes of diagrams are identified. The model 
potential is a Dirac-Hartree-Fock V^~^ potential, and B-spline basis functions in a cavity of finite 
radius are employed in the numerical calculations. The Breit interaction is taken into account 
through second order of perturbation theory and the lowest-order Lamb shift is also evaluated. 
Sample calculations are performed for berylliumlike ions with Z = 4-7, and for the magnesiumlike 
ion P IV. The third-order energies are in excellent agreement with measurement with an accuracy 
at 0.2% level for the cases considered. Comparisons are made with previous second-order MBPT 
results and with other calculations. The third-order energy correction is shown to be significant, 
improving second-order correlation energies by an order of magnitude. 
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I. INTRODUCTION 

The development of relativistic MBPT in recent 
decades has been motivated in part by the need for ac- 
curate theoretical amplitudes of parity non-conserving 
(PNC) transitions in heavy monovalent atoms such as 
cesium and francium. Applications of the theoretical 
methods developed to treat atomic PNC include support 
of atomic clock development, tests of QED in precision 
spectroscopy of highly-stripped ions, searches for time 
variation in the fine-structure constant, and providing 
precise astrophysical data. 

Although nonrelativistic studies QL lljl H H HQ 
and relativistic HF calculations [sL 13. Il0|| for divalent 
atoms and ions have been done for many years, only re- 
cently have relativistic many-body calculations been re- 
ported. As examples, we note that all-order relativistic 
MBPT calculations for transitions in berylliumlike ions 
with Z = 26 and 42 were carried out by Lindroth and 
Hvarfner while large-scale configuration-interaction 
(CI) calculations for transitions in C III were performed 



'Electronic address: hhol@nd.odu present address: NCTS 
(Physics Division), 101 Sec. 2 Kuang Fu Road, Hsinchu, Taiwan, 
R.O.C. 

tElectronic address: johnson@nd.eduj] 

Electronic address: steven.blundell @cea.trl 
§ Electronic address: msafrono@ udel.edul 



by Chen et al. J^]. Relativistic many-body calculations 
for magnesiumlike ions include the CI calculations of 
states in the n = 3 complex by Chen and Cheng [l3l |. 
and mixed CTMBPT calculations of excitation energies 
in neutral Mg by Savukov and Johnson [l^ . 

Second-order relativistic MBPT was applied to Be-like 
ions by Safronova et al. [l5l | and energies were found to 
be accurate at the 2% level. In this paper, we extend rel- 
ativistic MBPT for divalent atoms and ions to third or- 
der. We give a detailed treatment of the two-body terms 
here; the one-body terms are identical to those for mono- 
valent systems and are discussed in detail by Blundell 
et al. [ig . The long-range goal of the present research is 
to extend the relativistic singles-doubles coupled-cluster 
(SDCC) formalism to atoms and ions with two valence 
electrons. The present calculations permit us to identify 
and evaluate those third-order terms missing from the 
SDCC expansion. 

II. THEORETICAL METHOD 

The model potential for our MBPT calculation is 
the Dirac-Hartree-Fock V^^'^ potential. The Dirac- 
Coulomb-Breit Hamiltonian for an A'^-electron atom is 
H = Hq + V, where 

N 

Ho = ^[hu{i) + u(ri)], 

i=l 



2 



and the Dirac Hamiltonian is 

/iD = ca-p + f3c^ + Kuc(''), 

where Vnuc is obtained assuming a Fermi nuclear charge 
distribution. All equations are in atomic units. The per- 
turbation is 



V = 



^ / 1 



N 



configuration-weight vector. For simplicity, matrix ele- 
ments of the effective Hamiltonian are given for the un- 
coupled states, \vu)) = a^ai, |0) and \v'w'). In third or- 
der, each element consists of 12 one-body and 23 two- 
body terms. They represent a total of 84 one-body and 
218 two-body Goldstone diagrams. The multiplications 
of Clebsch-Gordan coefficients and the summations over 
magnetic quantum numbers are carried out during angu- 
lar decomposition. Only the two-body part of the third 
order is discussed here as the one-body part and complete 
second-order results are already presented in Refs. [l^ 



where bij is the Breit interaction 



1 

2^ 



In the no virtual-pair approximation, the excitations are 
limited to positive-energy eigenstates of H. 0, 0, 0, 

The eigenstates of a divalent system having angular 
momentum (J, M) are described by the coupled states 

^jm{vw) ^ ri^-ui ^ {jvjw,mvm^\JM) alal,\0) , (1) 



where |0) represents the ground state of the ionic core 
and 77t,u, is the normalization constant 



1 for V ^ w 
for V ^ w 

v2 



Here v and w specify the corresponding one-electron 
states with quantum numbers {ny,l^, j^,my) and 
{riw, lw,jw, rriio)- The model space P is defined by the set 
of total angular-momentum states the model-space 
projection operator is 

P= XI \^Jm{vw)) {^jMivw)\. 

JM 

(t,<™) 

The orthogonal-space operator Q is simply 1— P. 

The wave operator Q is found by solving the general- 
ized Bloch equation 

[fi, Ho] p = (vn - npvn)p. 

The effective Hamiltonian is given in terms of the wave 
operator 



H^g ^ PHaP + PVnP. 



A. Effective Hamiltonian 



We first find the configuration-weight vector by di- 
agonalizing the Hamiltonian H — H^J^ + us- 
ing total angular- momentum eigenstates 1^ as basis. 
Higher-order energies are obtained by operating the ef- 
fective Hamiltonian of the corresponding order on the 
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and [3. 



Z 

Si 

S2 



= E 



Do 



D9 



QcdvjvQahcdQru' v' ah 



= E 



QacmwQmhacQv' w' vb 



{£bv 



I {£a. 



X (^1 -\- \v -^-^ W , v' ^ w''^ ) + C.C. 

QahwmQmchvQv' w' ca 



E 



53 - E 



X (^1 -\- \v ^ w, v' ^ w'l ) + C.C. 

QahwmQw' cah9v' mvc 



5* - E 



X (1 — [u <-> to]) (1 — [w' It;'] ) + C.I 

QabwmQv' cvbQw' mac 



{£ac 



ahem 

x(l — [v ^ w\) (1 — \^v' ^ '^''\) 

7-1 \ ^ QahmnQmn'whQv' w' va 

= ~ 1 ri r 

ahmn 

X (1 + [li ^ U!, i;' -^-^ w''\ ) + C.C. 



= E 



' QahmnQmnvwQv' w' ah 



+ 



+ C.C. 



D3 = 
D4 = 

Ds = 
De = 



E' 



QahmnQw' nahQv' mvu 
(^^ahvw ^mnv' w'^ 



+ 



E 



QahmnQw' nwhQv' mva 



(1 + [f ^ v' ^ It;'] ) + C.C. 

(1 — [i; ^ w]) + C.C. 



E 



Qw' amnQnbawQv' mvb 
(^£avw ^mnv'') (,^bv ^-mu') 

x{l — [v ^ w]) (1 — [v' ^ It?'] ) + C.C. 

Qw' amnQnbvwQv' mab 

, \£avw ~ ^mnv' ) y^ah ^mv' ) 
abmn ^ 

X {1 + \v w , v' "r-^ w''\ ) + C.C. 

^ QabmwQw'mhnQv'nva 

, [£ahv ^mv' w' ) [^av ^nv' ) 
abmn 

X (1 — [i? ^ w\) (1 — \v' ^ w']^ ) + C.C. 

Qw' hwnQnahmQv' mva 



= E 



x{l — [v ^ w]) (1 — [v' ^ w']) 

QahwmQv' mvTiQw' nah 



E 



(^£abv £mv' w'^ (^£ab ^-jiw'^ 



x(l — [v ^ w]) (1 — [v' <^ w']) 



Ti 

Q 

B3 



E 



Qw'bmnQv' avbQmnwa 

[Sbv 



x(l — [f ^ w]) (1 — [u' ^ w']) 

' Qw' anmQmnarQv'rvw 



E 

x(l 

(^avw ^ S'mnt'') (^av ~ ^ru') 

amnr 

X (1 — [f ^ w]) (1 — [i;' ^ w''\ ) + C.C. 



X (^1 -\- [v ^ w, v' ^ w'] ) + C.C. 

Qw' amnQmnwrQv'rva 



- E' 



Qw' anrQv' rmaQmnvw 
(^Savw ^nrv'') {_^vw £-mn) 



amnr 

X (1 + [u ^ ItJ, t;' 1/;'] ) + C.I 
Qw' arnQv' nvmQrmwa 



E 

amnr 

x(l 

E 



x(l — [f <^ t/;]) (1 — \v' ^ w']) 

' Qv' w' rsQrsmnQmnvw 

Qw' amnQnmaxQv' xvw 



-E 

amn; 

E 



i^avw ^mnv'^ (^ax ^mji^ 

X (^1 -\- \v w , v' -^-^ w''^ ) 

Qv' axmQw' myaQxyvw 



{^avw ^mxw'^ {^a 



y '^mw' 



X (^1 -\- \v ^ w , v' -^-^ w''^ ) 

' Qv'w'mnQmnxyQxyvw 



E 



where giju = (ij \i'i2 | fc^) is the Coulomb matrix element 
defined in terms of single-electron functions as 



gijki 



ri - r2 



■0l(ri)0](r2)0fc(ri)0i(r2), 



and gijki = giju - gijik- The notation eiju = £i + £j + 
Ek + Si, etc. for the sum of single-electron eigenenergies 
has also been used. The third-order terms are arranged 
by the number of excited states in the sums over states. 
Zero, single, double, triple, quadruple excited-state terms 
are designated by the letters Z, S, D,T,Q. Terms associ- 
ated with backwards (folded) diagrams are designated 
by B. Backwards diagrams are unique for open-shell 
systems and exist only in the third or higher order of 
MBPT. The summation indices {a,b,c,d) refer to core 
states, (m, n, r, s) refer to excited states and, in back- 
ward diagrams, indices (x, y) refer to valence states. The 
primes above the summation signs indicate that excited 
states (m, n, r, s) are restricted to the orthogonal space 
Q only. This restriction applies only to the term with 
denominator l/(e«u, — £mn) in D2 and to the term with 
denominator l/{£vw—£mv') in D3. The c.c. denotes com- 
plex conjugate. The conjugate diagrams are obtained by 
a reflection through a horizontal axis, with the initial 
and final states switched (vw) <-> {v'w'). Direct dia- 
grams of the two-body terms are shown in Fig. ^ Tech- 
nically, there are subtle changes in energy denominators 



in going from the term presented to its c.c. counterpart. 
These changes can be deduced by re-drawing the diagram 
upside-down and reading off the new denominators. 

In Di and D4 (but not in D2 and D^) we combined 
diagrams associated with double excitations that have 
the same numerators but different denominators using 
the formula 

1 11 

{A + B)A ^ {A + B)B ^ AB' 

Diagrams Di - Di are special in the sense that two order- 
ings are possible. The ambiguous vertices are labeled by 
crosses in the diagrams. Many of the third-order terms in 
the two-body part have external exchanges and complex 
conjugates so each diagram illustrated in Fig. has from 
one to eight variants. The largest fraction of computer 
time is spent evaluating the term Q and most of the re- 
mainder is spent on repetitive evaluation of terms Df^ , Tf^ 
and their variants. Angular decompositions of the direct 
formulas are listed in the Appendix. 

III. APPLICATION AND DISCUSSION 

As a first illustration, we apply the theory described 
above to obtain energies of the ground state and ex- 
cited states in the n = 2 complex for Be-like ions. In 
Table HJ we give a detailed breakdown of the contribu- 
tions from first-, second- and third-order perturbation 
theory, together with corrections from the Breit interac- 
tion and from the Lamb-shift, for excitation energies of 
(2s2p) ■^-Po,i,2 states of Be-like ions. The experimental en- 
ergies are taken from the NIST database for atomic spec- 
troscopy |23|. Energies represent the lowest-order 
energies including the Breit correction. Lowest-order 
Lamb shifts -BLaiinb obtained following the method de- 
scribed in Ref. 23] . We find that the residual differences 
between calculated and measured energies AE decrease 
rapidly with increasing Z. This is expected since MBPT 
converges better for charged ions than for neutral atoms. 
In fact, for highly-charged ions, correlations are expected 
to decrease approximately as where n is the order 
of perturbation theory |2J|. On the other hand, QED 
effects (Lamb shifts) become more important along an 
isoelectronic sequence. The results in Table confirm 
both of these trends. Results of our calculations of ex- 
citation energies of all levels in the n = 2 complex for 
Be-like N (NIV) are presented in Table HH and are seen 
to be in agreement with measurement to parts in 10*. 

As a more involved example, we give a complete break- 
down of contributions to energies of low-lying states in 
the n ~ 3 complex for the Mg-like ion P IV in Table IIIII 
For both N IV and P IV, correlations are seen to account 
for about 10% of the total energies. For Be-like ions, the 
third-order correlations are only 7-10% of the second- 
order correlations. The second-order correlation energies 
are an order of magnitude smaller than the correspond- 
ing DHF energies. The Breit correction B'-^\ which is 
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TABLE I: Comparisons of third-order energies (cm^^) of the 
triplet (2s2p) states of Be-like ions Z—A-l with measure- 
ments are given, illustrating the rapid (l/Z^) decrease of the 
residual correlation corrections with increasing Z. A break- 
down of contributions to the energy from Coulomb and Breit 
correlation corrections and the Lamb shift is given. 
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23607.9 


39116.7 


54204.5 


69072.2 


£(2) 


-3114.3 


-2583.2 


-2344.0 


-2201.5 


5(2) 


-1.7 


-3.7 


-6.7 


-10.5 


£(3) 


473.5 


598.0 


412.9 


294.9 


^-Lamb 


-0.9 


-3.3 


-8.4 


-17.7 


Etot 


20964.6 


37124.4 


52258.4 


67137.3 


-S'expt 


21978.3 


37336.7 
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67209.2 


AE 


-1014 


-212 
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-2582.2 


-2342.2 
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-5.6 
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597.8 


412.7 


294.7 
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-8.3 


-17.5 
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67200.7 
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37342.4 
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67272.3 
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£(2) 
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-2338.4 


-2192.6 


5(2) 
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0.5 


0.4 
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£(3) 


473.2 


597.4 


412.3 


294.1 


-ELamb 


-0.8 


-3.2 


-8.2 


-17.1 


-Etot 


20968.0 


37147.3 


52339.1 


67345.2 


-Eexpt 


21981.3 


37358.3 


52447.1 


67416.3 


AE 


-1013 


-211 


-108 


-71 



obtained by linearizing the second-order matrix elements 
in Breit interaction, is also small for such lightly-charged 
ions. 

Our calculations for Be-like ions are able to produce 
results accurate to order of ten cm~^. This shows that 
the third-order energy correction is very important for 
divalent ions. By comparison, the second-order results 
of Safronova et al. [is'l for NIV agree with experiment 
at the level of a few hundred cm~^. It should be noted 
that results from the CI-I-MBPT method 14] mentioned 
in the introduction are consistently more accurate than 
the present third-order results. However, the CI-I-MBPT 
calculations contain a free parameter in the energy de- 
nominators that is adjusted to give optimized energies. 
In contrast, third-order MBPT is completely ah initio. 
For the C III ion, the large-scale CI calculations of Chen 
et al. |l2l | mentioned in the introduction also give tran- 



sition energies accurate to better than a hundred cm"-'^ 
on average. Those large-scale CI calculations are also ah 
initio and have about the same accuracy as third-order 
MBPT for states in C III. 

For the P IV ion, our results are in good agreement 
with experiment, the average dis crep ancy being several 
hundred cm~^. Chaudhuri et al. p5| employed an ejfec- 
tive valence shell Hamiltonian (EVSH) to calculate en- 
ergies of Mg-like ions and obtained results for P IV hav- 
ing a discrepancy of about a thousand cm^^, which is 
somewhat larger, but comparable, to the accuracy of our 
third-order calculations. 

It is informative to analyze our results in terms of di- 
agrams. The relative contributions of the third-order 
two-body terms for the ground-state energy of a typi- 
cal member in the Be sequence and the ion P IV are 
summarized in Tables IIVI and Dominant classes of 
diagrams are Q and -B3. This is understandable since 
Q are quadruple excited-state diagrams and involve no 
core excitations. The quadruple excitation diagrams are 
entirely due to valence-valence correlation effects; they 
are expected to be large because of the strong repulsion 
of the outer valence electrons. Class -B3 are backwards 
diagrams, which are characteristic of open-shell systems. 
As shown in Fig. ^ this class is also associated solely 
with valence-valence correlation. The two classes Q and 
-B3 tend to cancel each other as there is an extra phase of 
-1 associated with backwards diagrams. It is interesting 
to note that even after subtraction of the contributions 
from Q and .63, their difference is still larger than the 
contribution from any other single class of diagrams for 

cm. 



IV. CONCLUSIONS 

The accuracy of third-order MBPT results is at 0.2% 
level for lightly-charged ions of both Be and Mg isoelec- 
tronic sequences. This level of accuracy is comparable 
or superior to the two ah initio methods mentioned in 
Sec, mil A complete third-order calculation is important 
to understand the relative importance of different con- 
tributions to energies of divalent systems. The folded 
diagrams as well as the quadruple excitation diagrams 
are significant. The dominant role of these two classes 
of diagrams is attributed to the strong correlation of the 
two valence electrons. This conclusion is useful for work- 
ers developing combined CI-MBPT methods which in- 
clude dominant third-order diagrams. It is also helpful 
for researchers setting up SDCC calculations as they try 
to classify and account for the contributions from the 
third-order diagrams associated with omitted triple exci- 
tations. Although one might expect a complete fourth- 
order calculation for divalent systems to improve the ac- 
curacy of the present calculations still further, it is un- 
likely that such a complex calculation will be carried out 
in the near future. 
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TABLE II: Third-order energies (cm ^) of states in the n = 2 complex of the Be- like ion NIV, including corrections for the 
Breit interaction and the Lamb shift. 
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TABLE III: Comparison with measurement of theoretical en- 
ergies (cm^"*^) of some of the low- lying states in the n = 3 
complex of the Mg-like ion P IV, including a breakdown of 
contributions from Coulomb and Breit correlation corrections 
and the Lamb shift. 
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B(2) 


-0.9 


0.3 


1.3 


£;(3) 


807.4 


807.6 


807.6 


-E'Lamb 


-21.1 
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-20.5 


Etot 


67917.1 


68146.0 
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165654.0 


194591.8 


AE 


41 


3 


294 



TABLE V: Relative contributions of third-order two-body 
terms for P IV. 
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TABLE IV: Relative contributions of third-order two-body 
terms for C III. 
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APPENDIX 



Angular decompositions of direct formulas for the 
third-order two-body part are presented. 



^ XLi{cdwv)XL2{abcd)ZL3 (lo'v'ab) 

abed 

J jc jd [ \ ja jb J [ ) J ja jb 



Si = 



Ll jv jw J 1^ jd jc L2 j [ L3 j^i j„ 

1 y~^ XL{acmw)ZL{mbac)Zi^i (v'w'vb) 

abcm 



[L] y L' jb jv 
ZL{abwm)ZL{mcbv)Z]^i (v'w'ca) 



{Sac ^v'w'^ (^^abv ^mv'w'^ 



E 

LL' 
abcm 

1 \ J ja jc [ \ J jc ja 



m \ L jy jvu \ \ L' jw' jv' 

Xli {abwm)ZL2 (w'cab) Zl^ {v'mvc) 



\J+j^+ju, J -^1 -^2 L3 
Jc Jm Jb 



x(-l) 

jw jw' ^3 I J J jv' j 



L'z Ll ja j y Lz jw jv 

ZLiiabwm)ZL2 (v'cvb) Zls (w'mac) 



abcm 



x(-l) 



jv' jv 



L2 jw jv 
jw jw' L2 \ \ L\ L2 L3 



Di = 



L3 Ll ja ) { Jc Jm Jb 

1 ^ XL{abmn)ZL{rnnwb)Zj^i (v'w'va) 

\Jw^ ^ {^bw £mn) {Sav ^v'w'^ 

abmn 

^^_-|^^''+i'+Ja+Ji,+Jm+Jn+J,„/+J» 



x<5,„,„,-— 



J Jv' ju 



[L] L' 



Ja Jv 



Dn 



E' 

abmn 



XLi{abmn)XL2{T^nv'w)ZL3 {v'w'ah) 



(^£abvw 



■'mnv' w 



Dz = 



(^£vw Smn) i^Sab ^v'w'^ _ 

ja jb •^11 J jm 3n\\ J ja jb 

in jm Ll j [ L2 jw jv ] \ L3 jw' jv' 

1 ' Zi^(abmn)XL (w'nab) Zy {v'mvw) 



LL' 
abmn 



. i^ab ■ 



1 



') 



vw '-■mv 



)J 



x(-l) 



J+L'+ja +jb +im +3n +3v 



D4 



x5i 



J jv' jn 



[L] I L' j^ jv 
Zh(phrnn)Zh (w'nwb) Zl (v'mva) 



X ( — 1 ) ^ +-^^+Ja+J6 +Jm +Jn +J^(,/ +Jt> 

y J jw' jv' 

[L]'\L jv jw 

Zl {w'ainn) ZL{rihaw)Z (v'rnvb) 



LL' 
abmn 



(^Savw ^mnv'^ i^^bv ^'i 



') 



J jv' jw' 



De= y: 



jw' jw L 

[L] \ L' jw jv J \ jb jm L 

Zli {w'amn) Zl2 {nbvw)ZL3 {v'mab) 



(^Savw ^mnv'^ {^ab ^mv') 



^ |^_-|^'jl + i'l+-£'2 + i'3+J„'+Jt) 

J jn jb 

jw' Ll jm 
jv' ja L-i 

Zli {abmw)ZL2 {w'mbn) Zl^ (v'nva) 



J jn jb 
L2 jw jv 



- E 



L\L2L^ 
abmn 



{^abv ^mv'w'^ (^^av ^ 
jw jw' L3 



= E 



Ll L-z La 

Jn Ja Jm 

J jv' jw' 
L2 Ll jb J 1^ ^-3 jw jv 

Zl (w'hwn) ZL(nabrn)ZL (v'rnva) 



x(-l) 



[Lf 



i^Sbv ^nv'^ (^£av ^mv'') 



J+L-\-ja+3h-^jm+3n-\-3^t-\-3v 



J jw' jv' 
L jv jw 
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Ti 



T3 



'1^2 

abm 

x(-l) 



XLi{abwm)ZL2 {v'mvn) Zl^ {w'nab) 



^1 ^2^3 

abmn 



1) {£ab ■ 



') 



J+jv+jw J ^2 L:t 



L2 



J 3v' ju 



= - E 



L3 Ll ja ) { L2 jw jv 

Xli {w'bmn) (v'avb) ZL^imnwa) 



{Sbv 



'){ea 







abmn 
X(-1)'^+J"+^' 



Li L2 L3 

Ja Jn Jb 

jw' jw \ \ 'J jv jw' 

L3 Lj jm j y L2 jw jv 

1 ' Xl {w' anm) ZL{mnar)Zy (v'rvw) 

TiTl ^ 



x(-l) 



= E 



nnr 

J+l + L'+ja + Jm + Jn +jv 
J jv' jw' 

Xli {w'amn) ZL2i'mnwr)ZL;i {v'rva) 



x(-l) 



J+jv+jw I Li L2 La 

Jr Ja Jn 



J jv' Jw' 
L2 Ll jm J I ^3 

Zl (w'anr) Zl {v'rma) Z^r {mnvw) 



amnr 

^^_-|^^l+i''+ja +Jm + Jn + Jr + J^ / + j\; 



J j 



n Jm 



J j 



m Jn 



W I L j„/ j„/ J 1^ i' jw jv 

Zl^ {w'arn) Zl2 {v'nvm) ZL^irmwa) 



L-^ L'2 

amnr 



X{-1)' ' 

jw' jw L2 



J+Li+L2+L3+j^+j„ J Ll L2 L3 
Jm Ja Jn 



^1^2^3 
mnrs 

x(-l) 



J jv' jw' 
L3 Ll jr J ^2 jw jv 

f Xli (v'w'rs) XL2{^smn)ZL3{mnvw) 

J+Li+L2+L3+j„^+j„+jr+js+j^l+jv 



J jr js 
Ll jw' jv' 



jr js \ } j"'' j^ 

jn jm L2 I [ L3 ju, jv 



Bi = 



1 Z L {w' amn) X Linmax) Z 1^1 (v' xvw) 

x(-l)^ 



[jiti'] ^ {^avw ^mnv' 

amnx 



){ea 



J jv' ju 



[L\ 1 L' ju 



Jv 



B2 



Zl {v'axm) Zl {w'rnya) Zl' (xyvw) 



LL' 
amxy 



B3 = 



) 



m [ L j„/ ju' ^ [ L' jv, Jv 
y-\ ' Xli {v'w' mn) XL2(rnnxy)ZL3{xyvw) 



■^^_-^^i+L' +ja+jm+jx+jy+j„l+jv 
1 I J jx jy 11 J jx jy 



-mn) {,^xy 



mnxy 

J+Li+L2+L3+j^+j„+j^+jy+j^,+j^ 



X(-1) 

J jm jn 
Ll jw' jv' 



jm jn J 
jy jx L2 



J jx jy 
L3 jw jv 



The effective interaction strength 

XL{ijkl) = (-1)^ (i ||C^|| k) {j I) RLiijkl), 

is independent of magnetic quantum numbers. The re- 
duced matrix element of the tensor is 

{i\\c^k) = {-iy^+i.M\M(^^\ I ^^^u^k,ik,L), 

where [ji] = 2ji + 1 is the occupation number of shell i, 
and 



W {k,lk,L) 



1 \i li + Ik + L \s even 



\ili + lk + L \s odd 
The Slater integral RL{ijkl) is 



dridr 



^L+l 



[Pi{n)Pk{ri) + Q, {n)Qk (ri)] 

X [Piir2)Pi{r2) + Qj{r2)Qi (ra)] • 



The quantity ZL{ijkl) is defined by 



ZL{ijki) = XL{ijki) + [L]y2h^ ^! ^XxL'iijik). 
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